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ABSTRACT 



The transit spectrum of the exoplanet HD 189733b has recently been obtained between 0.55 and 1.05/xm. Here we present 
an analysis of this spectrum. We develop first-order equations to interpret absorption spectra. In the case of HD 189733b, 
we show that the observed slope of the absorption as a function of wavelength is characteristic of extinction proportional 
to the inverse of the fourth power of the wavelength (oc A~ 4 ). Assuming an extinction dominated by Rayleigh scattering, 
we derive an atmospheric temperature of 1340 ± 150 K. If molecular hydrogen is responsible for the Rayleigh scattering, 
the atmospheric pressure at the planetary characteristic radius of 0.1564 stellar radius must be 410±30mbar. However 
the preferred scenario is scattering by condensate particles. Using the Mie approximation, we find that the particles 
must have a low value for the imaginary part of the refraction index. We identify MgSiC>3 as a possible abundant 
condensate whose particle size must be between ~10 -2 and ^lO" 1 fim. For this condensate, assuming solar abundance, 
the pressure at 0.1564 stellar radius is found to be between a few microbars and few millibars, and the temperature is 
found to be in the range 1340-1540 K, and both depend on the particle size. 

Key words. Stars: planetary systems 
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1. Introduction 

Transiting extrasolar planets, like HD 189733b offer unique 
opportunities to scrutinize their atmospheric content (e.g., 
Charbonneau et al. 2002, Vidal-Madjar et al. 2003). In par- 
ticular, primary transits can reveal tenuous quantities of 
gas or dust by spectral absorption leading to variations in 
the apparent planet size as a function of the wavelength. 
HD 189733 b is presently the best target for such studies 
because of its nearby and bright host star, its large absorp- 
tion depth in its transit light curve (Bouchy et al. 2005), 
and a short period (Hebrard & Lecavelier des Etangs 2006). 

Using transit observations with the ACS camera of 
the Hubble Space Telescope, the apparent radius of 
HD 189733b has been measured from 0.55 to 1.05 microns, 
thus providing the "transmission spectrum" of the atmo- 
sphere (Pont et al. 2008). In this wavelength range, strong 
features of abundant atomic species and molecules were pre- 
dicted but not detected, leading Pont et al. (2008) to the 
conclusion that a haze of sub- micron particles is present in 
the upper atmosphere of the planet. 

Here we use these measurements of the planet radius 
as a function of the wavelength from 550 to 1050 nm. In 
Sect. [5] we describe how a transmission spectrum can be in- 
terpreted to derive basic quantities. We then estimate the 
temperature of HD 189733b from fit to the ACS spectrum 
and show that the observed absorption can be caused by 
Rayleigh scattering (Sect. [3]). Several species that are pos- 
sibly responsible for the Rayleigh scattering are discussed, 
and the corresponding pressure at the absorption level are 
estimated in Sect. |U 



2. Interpreting a transit spectrum 

2.1. Planet radius as a function of atmospheric content 
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Fig. 1. Plot of r eq as a function of the ratio of the planet 
radius to the atmosphere scale height. 



To interpret transit spectra, one can use detailed atmo- 
sphere models that include temperature, pressure, and com- 
position as a function of the altitude and numerically inte- 
grate the radiative transfer equations to obtain a theoretical 
spectrum to compare to the observations (e.g., Ehrenreich 
et al. 2006). Here we propose another approach by deriving 
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first-order equations to obtain a better feeling for the basic 
quantities that can be obtained from these measurements. 

Following Fortney's (2005) derivation of the path 
length, the optical depth, r, in a line of sight graz- 
ing the planetary limb at an altitude z is given by 
t(X,z) « a(X)n(z)y/ 2?r R p i anet if , where R p i an et is the as- 
sumed planet radius, if the atmosphere scale height, and 
n(z) = n( z= Qjexp(— z/H) the volume density at the alti- 
tude z of the main absorbent with a cross section cr(A). 
For a temperature T, if is given by if = kT/fig, where 
/i is the mean mass of atmospheric particles taken to be 
2.3 times the mass of the proton, and g the gravity. Using 
planetary transit measurements, a fit to the light curve pro- 
vides the ratio of the effective planetary radius as a func- 
tion of wavelength to the stellar radius: R p (X)/R*. Here 
we define r eq by the optical depth at altitude z eq such 
that a sharp occulting disk of radius R p i an et + z eq pro- 
duces the same absorption depth as the planet with its 
translucent atmosphere; in other words, r eq is defined by 
R P (X) — Rpianet + z[r = T eq ) . Using a model atmosphere, 
and numerically integrating over the whole translucent at- 
mosphere, we can calculate the effective planet radius, and 
then obtain the corresponding optical thickness at the effec- 
tive radius, r eq . For various atmospheric scale heights, we 
calculate r eq by numerical integration. For a wide range of 
atmospheric scale height, provided that R p /H is between 
~30 and ~3000, the resulting r eq is roughly constant at 
a value r eq w 0.56 (Fig. HJ. In the case of HD 189733b, 
Rp/H varies between 280 and 560 when the temperature 
varies from 1000 to 2000 K; therefore, the approximation 
of a constant r eq at 0.56 fully applies. This demonstrates 
that, for a given atmospheric structure and composition, 
estimating the altitude at which t = r eq = 0.56 is all that 
is needed to calculate the effective radius of the planet at 
a given wavelength. 

2.2. Consequences 

For a given atmospheric structure and composition, the ef- 
fective altitude of the atmosphere at a wavelength A is cal- 
culated by finding z(t — T eq , A), which solves the equation 
t(A, z) — T eq . Using the quantities defined above, the effec- 
tive altitude z is given by 



(A) = if In ( £,absPz=o<7abs(X)/T eq x J 2n R p I kT \xg I , 



where a aba and £ Q bs are the cross section and abundance of 
the dominant absorbing species. 

If the variation of the cross section as a function of wave- 
length is known, the observation of the altitude as a func- 
tion of wavelength allows the derivation of if and, therefore, 
of the temperature T, given by: 



T = — ( dlna \ 1 dz ^ 
k 1 dX J dX ' 



(2) 



Using Eq. [TJ the partial pressure of the main absorbent 
at the reference altitude is estimated by 



£,ab s Pz=o = T eq /<j abs (X z = ) x yj kTfig/2nR p , (3) 

where X z= o is the wavelength at which the effective plan- 
etary radius corresponds to (or is used to define) z = Q. 



We note that there is a degeneracy between the abundance 
and the total pressure in the atmosphere. From the mea- 
surement of the effective radius using transit spectroscopy, 
one needs to assume an abundance of the absorbent to de- 
rive the pressure, or alternatively to assume a pressure to 
derive the abundance. 

In summary, the effective planetary radius is character- 
istic of the pressure and abundance, and the variation in 
this radius as a function of wavelength is characteristic of 
the temperature. 



3. The atmospheric temperature in HD 189733b 

In HD 189733b, the plot of altitude as a function of the 
wavelength shows an increase in absorption toward shorter 
wavelengths (Fig. [2]). Using Eq. [2] and assuming a scaling 
law for the cross section in the form a = ero(A/Ao)", the 
slope of the planet radius as a function of the wavelength is 
given by dR p /d In A = dz/d In A = aH. Therefore, we have 



aT = 



k d\n X 



(4) 



For HD 189733b, the slope is measured to be 
472km±42(stat.)±106(syst.) from 600 to 1000 nm (Pont 
et al. 2008). This corresponds to dR p /d In A w -920 km. 
We therefore obtain aT«- 5840 K. The temperature in 
HD 189733b is determined to be in the range 900-1500 K 
(Doming et al. 2006; Knutson et al. 2007); therefore, a is 
found to be close to a « —4, which is typical of Rayleigh 
scattering. 

Assuming a = —4 as expected for Rayleigh scattering 
and using Eq. [21 we derive a temperature 

T = 1460 ± 130(stat.) ± 330(syst.) K, 

where (stat.) and (syst.) correspond to statistical and sys- 
tematics error bars on dR p , as quoted in Pont et al. (2008). 

For condensates, the scale height if c is often found to 
be significantly smaller than the gaseous scale height (if c ~ 
if/3) (Ackerman & Marley 2001; Fortney 2005). If this were 
the case, because of the observed slope in the spectrum, the 
temperature (or a) should be about three times higher than 
found above, which seems unlikely. We therefore conclude 
that, if produced by dust condensates, the observed transit 
spectrum of HD 189733 b shows that these condensates must 
be unusually well- mixed vertically with the atmospheric gas 
with a similar scale height. 



4. Rayleigh scattering 

4.1. Molecular hydrogen 

One possible carrier of the Rayleigh scattering is the 
most abundant molecule, molecular hydrogen. The effec- 
tive altitude at which the Rayleigh scattering of molecu- 
lar hydrogen dominates only depends on the mean den- 
sity. Therefore, measurements of the altitude in the regime 
where Rayleigh scattering by H2 dominates over other ab- 
sorbents allows the determination of not only the tem- 
perature, but because the abundance of H2 is close to 1 
(?»2 ~ 1) a l so the total density and consequently the total 
pressure at the reference zero altitude. 
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Fig. 2. Plot of the altitude in the atmosphere of 
HD 189733 b corresponding to the measured planetary ra- 
dius as a function of wavelength. Here the zero altitude 
is defined for a planetary radius of 0.1564 times the stel- 
lar radius. The measured values with error bars are those 
of Pont et al. (2008). The fit to the data (thick line) is 
obtained assuming an extinction proportional to A~ 4 and 
for an atmospheric temperature of 1340 K. Fits using H2 
Rayleigh scattering or MgSiC>3 grains with sizes between 
0.01 and 0.1 /im give the same plots. Larger grain diame- 
ters are not consistent with the observations (e.g., 0.4/im, 
shown with dotted line). 

Following Eq. [3] the pressure Po at the altitude corre- 
sponding to the radius at wavelength Aq is 



Po 



eq/o-Q x ^JkTng/2nR p , 



where <tq is the Rayleigh scattering cross section at Ao- 
Using the refractive index of molecular hydrogen (r — 1) = 
1.32 x 10~ 4 , one derives a = 2.52 x 10 _28 cm 2 at A = 
750 nm (Bates 1984; Naus & Ubachs 2000). Using these 
numbers one finds Po ~ 400mbar. 

A fit with two free parameters (temperature and pres- 
sure) and assuming that Rayleigh scattering by molecu- 
lar hydrogen dominates leads to a very good fit to the 
data published by Pont et al. (2008) with a x 2 equals 9.8 
for 8 degrees of freedom. We thus obtain a temperature 
of T=1340±150K and a pressure at z — (defined by 
P p /P* = 0.1564) of P=410±30mbar. 

4.2. Condensate 

Another possible carrier of the Rayleigh scattering is haze 
condensate typically smaller than the wavelength. 

Assuming that the absorption measured between 550 
and 1050 nm is due to particle condensate, we calculate 
the transmission spectrum using the Mie approximation. 
In transmission spectroscopy, given the small solid angle 
defined by the star seen from the planet, we can use the 
single scattering approximation. The glow produced by off- 
axis scattering is also negligible (Hubbard et al. 2001). In 
this configuration, the Mie extinction efficiency (Q ex t) is 
used and calculated as the summation of the absorption 
efficiency (Q a bs) and scattering efficiency (Q sca )- 

Rayleigh scattering in the form A~ 4 is obtained 
only for particle sizes (a) much smaller than the wave- 
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Fig. 3. Plot of the temperature and pressure at zero al- 
titude as a function of the size of MgSiC>3 grains needed 
to obtain a satisfactory fit to the data as in Fig. [5] Here 
we assume solar abundance for magnesium atoms to derive 
the grain abundance. For particles of 0.03 fim, temperature 
is found to be the same as with H2 Rayleigh scattering : 
T=1340± 150 K. At 0.01 and 0.1 /im temperature increases 
slightly to T=1540±170K and T=1390±150K, respec- 
tively. For particles smaller than ~0.01 /im and larger than 
~0.1/xm, Rayleigh scattering does not dominate the ex- 
tinction (— a^;4), and the temperature must be higher to 
fit the slope of the measured radius as a function of wave- 
length. At 0.3 \i the fit is bad with % 2 =22.3 for 8 degrees of 
freedom; \ 2 increases steeply for larger sizes. 

length (27ra<cA). In this last case, the absorption and 
scattering efficiencies are given by the approximation: 
Q abs = (87ra/A)3((n 2 - l)/(n 2 + 2)), and Q sca = 
8/3(27ra/A) 4 3fi(((n 2 - l)/(n 2 + 2)) 2 ), where n = n + ik 
is the complex refraction index of the condensate. 

We see from the two last equations that Rayleigh scat- 
tering dominates the transit spectrum only if Q a bs *C Qsca- 
Assuming that the imaginary part of the refraction index is 
much smaller than the real part, we find that this condition 
occurs for particle size, a, larger than a minimum size: 



a > a„ 



= 0.331* 



nk 



(n 2 



1) 



For many plausible dust condensates (like Mg2SiC>4, 
MgFeSiC>4, O-deficient silicates, etc.), the imaginary part 
of the refraction index is about k ~ 0.1, and the real part 
is in the range n ~ 1.5 — 2.0 (Dorschner et al. 1995). Using 
these values, we find that the minimum particle size is in 
the range a m j„ = A/11 - A/7. This provides very few pos- 
sibilities for the particle size, which must be much smaller 
than A/27T and much larger than ~A/10. We therefore con- 
sider unlikely that the transit spectrum of HD 189733b can 
be due to particles with k ~ 0.1, which would require fine- 
tuning of the particle size for the Rayleigh scattering dom- 
inating between 550 and 1050 nm. 

A broader range of particle size is allowed if the imag- 
inary part of the refraction index is significantly smaller 
than ~0.1. In other words, the condensate must be trans- 
parent enough for the absorption to be negligible relative 
to scattering. Among the possible abundant condensates, 
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MgSiC>3 presents such a property (Dorschner et al. 1995; 
Fortney 2005). If MgSiC>3 contains no Fe atoms, the imagi- 
nary part of the refraction index is found to be k ~ 2 x 10~ 5 
in the considered wavelength range. Using this value, we 
find a m i n /X = 8 x 10~ 3 . Therefore, we conclude that 
MgSiC>3 is the candidate species possibly responsible for 
the observed Rayleigh scattering in the HD 189733b spec- 
trum, and its size must be in the range 10 _2 -10 _1 microns 
(see upper panel of Fig. [3]). 

Assuming an abundance for the MgSiC>3 grains and 
using Eq. [3l we can evaluate the pressure at the effec- 
tive planetary radius if the Rayleigh scattering by these 
grains dominates the transit spectrum. For MgSiC>3, the 
refraction index is taken to be n~1.6 in the range 550- 
1050 nm (Dorschner et al. 1995). This gives a cross section 
a w 1528a 6 /A 4 . The grain abundance is obtained using so- 
lar abundance and assuming that it is limited by the avail- 
able number of magnesium atoms. Thus the grains have 
an abundance £ grai n = 2£ M gMMgSi0 3 /(/ 3 grain47r/3a 3 iVm p ), 
where ^Mg ~ 4 x 10~ 5 is the magnesium abundance, 
/iM g Si0 3 is about 100.4 times the mass of the proton (m p ), 
Pgrain = 3.2 g cm -3 is the grain density, a the particle size, 
and N Avogadro's number. Therefore we have £ gra in = 
1.0 x 10~ 15 (a/l^m)- 3 . Finally, at z = defined by the 
effective radius at A rj 750 nm, with T=1340K we obtain 
a pressure P = 2 x 10 _3 bar if the size of MgSiOs grains 
is about 0.01/im, and P = 2 x 10~ 6 bar if the size is about 
0.1/im. These values obtained with Eq. [3] are similar to 
those obtained with a general fit to the data (Fig. [3]) . It is 
noteworthy that these values are consistent with the con- 
densation pressure of MgSiC>3 which is ~10 -4 bar at 1400 K 
(Lodders 1999). 

Finally, Rayleigh scattering by MgSiOs is preferred to 
scattering by H2 to explain the observed spectrum. Indeed, 
molecular hydrogen requires higher pressure, at which the 
signature of abundant species like Nai, Kl, H2O, and possi- 
bly TiO should overcome the Rayleigh-scattering signature, 
except when they have anomalously low abundances. 

Because the cross section varies as oc a 6 in the Rayleigh 
regime, the scattering is largely dominated by the largest 
particles in the size distribution. Therefore, the typical sizes 
quoted above must be considered as the size of the largest 
particles in the distribution. 

5. Discussion 

This work assumes that the whole transit spectrum from 
550 to 1050 nm can be interpreted by the absorption by 
a unique species. Alternatively, there could be a combina- 
tion of absorptions by various species (e.g., Nai, Kl, and 
H2O) in different parts of the spectrum to mimic a spec- 
trum of Rayleigh scattering with the expected atmospheric 
temperature. This would be coincidental, but it cannot be 
excluded from the present data. This question needs higher 
resolution data to be solved. 

The temperature and pressure found above are mean 
values integrated along the planetary limb, where temper- 
ature is expected to vary from the pole to the equator and 
from one side to the other. A numerical integration of the 
absorption through the planetary atmosphere assuming a 
variation in temperature as a function of altitude of -60 K 
per scale height (Burrows et al. 2007) or as a function of 
longitude of 350 K over 120 degrees (Knutson et al. 2007) 



shows that temperatures found by solving Eq.[T]are slightly 
overestimated by no more than about 3% and 2%, respec- 
tively. 

If the measurements of planetary radius in the 550- 
1050 nm wavelength band are extrapolated to the infrared, 
where Spitzer measurements are available, we can ob- 
tain a minimum effective radius due to Rayleigh scatter- 
ing. From i? p /i?*=0.1564 at 750 nm, we extrapolate that 
Rp/R* must be larger than 0.1539, 0.1531, and 0.1526 
at 3.6 fim, 5.8 ftm, and 8.0 ^m, respectively (with ±0.0003 
(stat.) 0.0006 (syst.) error bars). These values are consistent 
with the Spitzer/IRAC measurements given by Knutson et 
al. (2007) and Ehrenreich et al. (2007), and with measure- 
ments by Beaulieu et al. (2007) at the limit of the upper 
error bars. 

Rayleigh scattering has also been proposed to explain 
the detection of polarized scattered light (Berdyugina et 
al. 2007). However, this detection concerns the high upper 
atmosphere and cannot be directly compared to the present 
work on scattering deeper in the atmosphere. 

Ground-based, high-resolution spectra allowed the de- 
tection of sodium in the atmosphere of HD 18973b (Redficld 
et al. 2008). Assuming a pressure at a given altitude, 
the measured absorption allows the determination of the 
sodium abundance. For that purpose, a fit of the sodium 
line shape is needed, and still has to be developed. 
Alternatively, by assuming the sodium abundance, mea- 
suring sodium absorption should constrain the pressure 
and allow for differentiating between the several carriers 
of the Rayleigh scattering proposed in the present paper. 
Presently, extinction high in the atmosphere by MgSiOs 
with its low condensation pressure of ^10~ 4 bar appears a 
natural explanation for the absence of broad spectral sig- 
nature and is the preferred scenario. 
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